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Abstract 

We study the electrophoretic separation of polyelectrolytes of varying lengths by means of end- 
labeled free-solution electrophoresis (ELFSE). A coarse-grained molecular dynamics simulation 
model, using full electrostatic interactions and a mesoscopic Lattice Boltzmann fluid to account 
for hydrodynamic interactions, is used to characterize the drag coefficients of different label types: 
linear and branched polymeric labels, as well as transiently bound micelles. 

It is specifically shown that the label's drag coefficient is determined by its hydrodynamic 
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size, and that the drag per label monomer is largest for linear labels. However, the addition of 
side chains to a linear label offers the possibility to increase the hydrodynamic size, and therefore 
the label efficiency, without having to increase the linear length of the label, thereby simplifying 
synthesis. The third class of labels investigated, transiently bound micelles, seems very promising 
for the usage in ELFSE, as they provide a significant higher hydrodynamic drag than the other 
label types. 

The results are compared to theoretical predictions, and we investigate how the efficiency of 
the ELFSE method can be improved by using smartly designed drag-tags. 



Introduction 

As known from experiments and theory, the free-solution mobility of a flexible polyelectrolyte 
chain does not depend on the chain length N (number of monomers) anymore if the chain is longer 
than a certain length A/po- The regime where > Npd is called free-draining regime. In this regime, 
the counterions influence the inter-monomer hydrodynamic interactions and allow the fluid to drain 
through the polyelectrolyte coil. The effective friction Feff becomes linear in the chain length, as 
does the effective charge QgfF for longer chains, which leads to a constant, length-independent 
mobility 

(1) 

i eff 

It was shown that attaching a suitable uncharged molecule to an electrophoresis target can re- 
store the size-dependent mobility and overcome the free-draining property of long polyelectrolyte 
chains. This method, which is known as end-labeled free-solution electrophoresis (ELFSE), 
is based on the alteration of the charge-to-friction ratio of the polyelectrolyte molecules by an 
uncharged drag label. 

The effect of the label can be compared to that of a parachute attached to a moving object. 
The additional friction provided by the parachute slows the object down. This effect is stronger 
for smaller molecules with a lower effective friction and smaller charge to pull the drag-tag, as the 



2 



Kai Grass et al. 



Optimizing ELFSE 



ratio between charge and friction is changed more drastically than for larger molecules. 

Since the method's introduction, finding suitable labels that provide a high hydrodynamic drag 
has been a major concern in this field. A larger hydrodynamic drag enables the separation of 
longer chain fragments, as the length-dependence of the electrophoretic mobility decreases with 
increasing polyelectrolyte length. When the additional friction provided by the drag-tag becomes 
negligible against the intrinsic effective friction of the polyelectrolyte, the chain becomes essen- 
tially free-draining again. Experimentally relevant is that the mobility of long polyelectrolyte 
chains should differ by a factor large enough to allow for accurately separating them, although 
their lengths only vary by a single monomer. The maximum chain length resolvable in this way is 
called the read length. 

In general, the drag labels can be chosen from a wide range of molecules but they have to fulfill 
certain requirements, such as being water-soluble at experimental conditions, having a unique 
attachment mechanism to the polyelectrolyte and showing minimal interaction with the walls of the 
capillary. The read length is optimised by choosing a large molecule that imposes a high frictional 
drag. However, to fulfill resolution requirements, the labels must remain perfectly monodisperse 
since polydispersity will effectively be like an additional source of diffusion that would broaden 
the peaks. 

As it poses an experimental challenge to produce large, monodisperse linear polymer labels, 
two recently proposed alternatives seem promising. Haynes et. al. proposed to use branched 
polymers with well-defined architectures. A first theoretical study on this method verified the 
approach and concluded that, even though a branched polymer is more compact and thus provides 
a smaller hydrodynamic friction for a given molecular weight than a linear polymer, this draw- 
back is offset by the monodispersity of the branched labels created by assembling shorter linear 
chains. Grosser et. al. ^^"^^ introduced nonionic surfactant micelles as drag labels with very large 
hydrodynamic friction. The inherent polydispersity of the micelles is overcome by using a PNA 
amphiphile that only provides a transient binding between the DNA fragments to be separated and 
the micelles. Each fragment attaches to a different micelle every couple of seconds, which results 
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in an averaging procedure over the course of the elution time that remedies the need for perfect 
monodispersity. 

The theoretical description of the methods discussed above is not complete and furthermore 
predicts behaviours that are difficult to test experimentally, such as end of chain effects, the hydro- 
dynamic deformation of the label in high fields, or the steric segregation between the label and the 
chain. Since it is not possible to visualize DNA-label molecules in the lab, computer simulations 
can support the understanding of the real physics as long as they include hydrodynamic interactions 
between polyelectrolyte, label and solvent, as well as account for the influence of the electrostatic 
interaction between the polyelectrolytes and its surrounding counterions. Of course, it is beyond 
the scope of this article to cover all previous predictions. 12,14,18-21 jjowever, we will demonstrate 
that it is possible to study these factors, and that the standard theory appears to be sufficient for the 
cases treated here. More cases will be studied in future papers. 

Since the ELFSE method overcomes the main drawback of ordinary gel electrophoresis, the 
long separation time due to the slow down by the applied gel matrix, it is a promising method on 
the way to faster sequencing methods and, as such, of especial interest to the community. 

In this paper we will use coarse-grained MD simulations to study the electrophoretic separa- 
tion of fully flexible polyelectrolytes of varying lengths by end-labeling. After introducing the 
simulation model, we confirm that the free-draining behaviour is correctly reproduced and test 
the standard theory for ELFSE by attaching an uncharged linear label. In Section , we introduce 
branches to the drag label and test the predictions made by. ^^'^^ From the branched label, we will 
go to micellar targets (Section ) and analyze the method proposed by. ^^'^^ We establish a relation 
between the average size of the micelle and its drag value. Our concluding remarks point out the 
efficiency of the ELFSE method and show the benefit of the different labeling methods. 
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Theory 

The theory for end-labeled free-solution electrophoresis is based on the interplay between hydro- 
dynamic and electrostatic forces, and it takes into account the stress that builds along the chain 
backbone. In general, it is assumed that the electrostatic and factional forces do not deform the 
hybrid molecule's random coil conformation nor its cloud of counterions. Because of these as- 
sumptions, the theory used here is valid for low velocities and weak electric fields. 

Neglecting molecular end-effects, the electrophoretic mobility jU = v/fi of the poly electrolyte 
with an attached linear drag-tag can be described in terms of the effective friction of the polyelec- 
trolyte Fpe, its effective charge Qpe and the hydrodynamic friction of the attached label Fl: 

2pe 1 

where jUq is the length independent free solution mobility without drag-tag. 

Equation 2 shows the importance of the ratio between Fpe and Fl. The electrophoretic mobility 
jU is a function of for a fixed Fl as long as FpE changes with and the ratio between Fl and 
FpE remains non-negligible. 

Since the electrophoretic friction coefficient FpE grows linear with the length of the polyelec- 
trolyte for long chains, as shown in a previous publication,^^ Equation 2 can be reformulated as 
follows: 



with a constant drag coefficient 



Here the ratio Fpe/N is the friction per monomer of the poly electrolyte. The chemistry and tem- 
perature dependent Ol is a measure for the difference in hydrodynamic properties between the 
polyelectrolyte and the label and represents the number of polyelectrolyte monomers that provide 
a hydrodynamic friction equal to that of the label. ^'^^'^^'^^'^^ 
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In order to characterize the effectiveness of an arbitrary (not necessarily Unear) label for size- 
separation. Equation 3 has been used to define this specific label property from the measured 
mobilities: 

= (5) 
where ttL is conveniently determined as the slope when plotting jUq/m versus l/NpE- 

^ = l + aL/N. (6) 



Simulation model 

We employ coarse-grained molecular dynamics (MD) simulations using the ESPResSo package 
to study the electrophoretic separation of fully charged linear flexible polyelectrolytes by end- 
labeled free-solution electrophoresis. The polyelectrolytes are modelled by a totally flexible bead- 
spring model as a set of spheres that represent the N individual monomers which are connected to 
each other by finitely extensible nonlinear elastic (FENE) bonds: 

with stiffness k = SOEq, and maximum extension R = 1.5c7o, and r being the distance between 
the interacting monomers. Additionally, a truncated Lennard- Jones or WCA potential is used for 
excluded volume interactions:^^ 

with the cut-off being rgut = V^cr, at which Uuir) = 0. 

The parameter Eq and Oq define the energy and length scale of the simulations. We use £0 = 
k^T, i.e. the energy of the system is expressed in terms of the thermal energy. The length scale 
(To defines the size of the monomers. We set <7o = 4.7A, which with the average bond length 
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along the PE chain being 0.91(7o represents a Unear monomer distance of approximately 4.3 A, 
the spacing of single-stranded DNA. Different polyelectrolytes can be mapped by changing Gq. 
Unless mentioned otherwise, all observables are expressed in terms of the simulation units Gq and 
eo, which we will not use explicitly from now on. 

Besides the dissociated counterions the system also contains additional monovalent salt. The 
counterions and the salt ions are modelled as charged spheres using the same WCA potential giving 
all particles in the system the same size. 

All chain monomers carry a negative electric charge q = —leo, where eo is the elementary 
charge. For charge neutrality, N monovalent counterions of charge +leo are added. Additional 
monovalent salt is added to the simulation. Full electrostatic interactions are calculated with the 
P3M algorithm using the implementation of Reference .^^ The Bjerrum length Ib = e^/Anek-QT = 
1.5 in simulation units corresponds to 7.1 A, the Bjerrum length in water at room temperature. 
This means that the effect of the surrounding water is modelled implicitly by simply employing 
the dielectric properties of water, using a relative dielectric constant of fa 80. The applied 
external field £ = 0.1 is represented by a constant force, acting on all charged particles. The field 
strength is chosen to be sufficiently small in order to not influence chain conformations or couterion 
distributions. 

The coarse-grained molecular dynamics model used in this article is extended by the inclusion 
of the three drag-tag labels investigated, as shown in Figure Figure 1 . The linear labels use the same 
flexible bead- spring model as the polyelectrolytes, but are uncharged. For the study of branched 
labels, flexible side chains of well-defined length are added to the linear label. The third kind of 
label, micellar drag-tags, is represented by a sphere of given radius whose surface is modelled by 
many small WCA spheres that are connected with each other by a network of FENE springs. The 
number of small spheres is defined by the radius of the large sphere to be modelled. This model 
has been successfully used to study colloidal electrophoresis. 

We include hydrodynamics using a Lattice Boltzmann (LB) algorithm that is frictionally 
coupled to the MD simulations via an algorithm introduced by Ahlrichs et al. The mesoscopic 
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Figure 1: (a) Polyelectrolyte with surrounding counter- and co-ions, (b) with linear drag-tag, (c) 
with branched polymeric drag-tag, and (d) with micellar drag-tag. 
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LB fluid is described by a velocity field generated by discrete momentum distributions on a spatial 
grid, rather than explicit fluid particles. We use an implementation of the D3Q19 model with a 
kinematic viscosity v = 1.0, a fluid density p = 1.0.^^ The resulting fluid has a dynamic viscosity 
T] = pv = 1.0. The space is discretized by a grid with spacing a = 1.0. The fluid is coupled to the 
particles by a frictional coupling with bare friction parameter Fbare = 20.0. Random fluctuations 
for particles and fluid act as a thermostat. The interaction between particles and fluid conserve 
total momentum and are proved to yield correct long-range hydrodynamic interaction between 
individual particles. 

The simulations are carried out under periodic boundary conditions in a rectangular simula- 
tion box. We investigate the behaviour of polyelectrolyte chains varying from AT = 20 to = 60 
monomers. The size L of the box is varied to realize a constant monomer density of npg — 10~^ 
which corresponds to a concentration cpE = 16mM. The same concentration is used for the addi- 
tional salt, resulting in a Debye length of Ad ~ 4.2. A MD time step Tmd = 0.01 and a LB time 
step TLB = 0.05 are used. After equilibration of 10^ steps, 10^ steps are used for generating the 
data. The time-series are analyzed using auto-correlation functions to estimate the statistical errors 
as detailed in Reference. Error bars of the order of the symbol size or smaller are omitted in the 
figures. 

The electrophoretic mobility is obtained by applying a constant electric field of reduced field 
strength £ = 0.1 that acts on all charged particles. The mobility is then given by direct measure- 
ment of the center of mass velocity v of the chain: 

" = T 

Before applying this method, it was ensured that the applied electric field strength E is small 
enough not to distort chain conformations or counterion distributions. Therefore, the system is in 
the linear response regime, i.e., the measured mobility does not depend on the magnitude of the 
electric field. 
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Up to ten independent simulations are carried out for each data point, taking between one day 
and two weeks on a single standard CPU^ depending on the chain length N and the t5^e of label 
investigated. 

Linear drag tags 

In this section, the simulation model is applied to the electrophoresis of polyelectrolyte chains 
with an attached linear polymeric drag-tag. The electrophoretic mobility for polyelectrolyte chains 
is determined with and without different labels, and the results are compared to the theoretical 
predictions. We also examine how the effective friction of the drag-tag is influenced by the intrinsic 
stiffness and the salt concentration in the solution. 

Testing the standard ELFSE theory 

First, the free-solution electrophoretic mobility without an attached drag-tag, jio, is determined, as 
shown in Figure Figure 2. The measured mobility does not depend on the chain length, as expected 
for longer free-draining polyelectrolyte chains. The average mobiUty is determined to be 

Aio = 0.147 ±0.002. (8) 

Additionally, in Figure Figure 2, the mobilities with attached drag-tags ranging from L = 30 
to L = 70 monomers are measured, and it is confirmed that a length-dependence is achieved and 
that the difference in mobilities, i.e., the selectivity of the separation, is better the longer the at- 
tached label is. Equation 6 is used to calculate the hydrodynamic drag coefficients as shown in 
Figure Figure 3, resulting in values from = 13.5 ± 0.4 for L = 30 to ttL = 25.7 ± 0.6 for L = 70. 

In the following, an expression for Ol based on the hydrodynamic size and shape of the label 
is developed. The hydrodynamic friction Fl of the uncharged label is related to the hydrodynamic 

^Dual Core AMD Opteron(tm) Processor 270 
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Figure 2: The free-solution electrophoretic mobility without label (black circles) shows no depen- 
dence on the chain length A'^. The free-draining mobility is jiQ = 0.147 ±0.002. The attachment 
of linear drag-tags to the end of the polyelectrolyte chains reduces the mobility and restore a A^- 
dependent behaviour. The label length L is varied from 30 to 70 monomers, with the largest label 
resulting in the strongest slowdown. 
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Figure 3: The hydrodynamic drag coefficient (Zl is given by the slope of the curve. For the linear 
labels, ttL ranges from 13.5 ±0.4 to 25.7 ±0.6. 
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radius by means of the Stokes relation: 

Tl = 6;r77i?h,L- (9) 

As in the previous chapter, the effective electrophoretic friction of the polyelectrolyte is expressed 
in terms of the free-solution mobility jUq and the effective charge Q^s: 

TPE = eeff/A^o- (10) 

Using the Manning prediction Qeff = (l/i^) for the effective charge, where the Manning con- 
densation parameter ^ = l^/b is a measure for the strength of the electrostatic potential of the 
polyelectrolyte, finally yields 

aL = jUo'^6;rT]i?h,L- (11) 
With the system parameter used here, ^ = 1.63, one obtains 

aL = (4.5±0.1)i?h.L- (12) 

Equation 1 1 will be shown to be valid for linear labels whose size is not exceeding the Debye 
length Ad. When the label size becomes larger, the friction of the label is not anymore directly 
related to the hydrodynamic radius, as the salt ions that penetrate the polymer coil influence the 
intermonomer hydrodynamic interactions and limit them to the electrostatic screening length. As 
for the polyelectrolyte itself, this screening length is of the order of the Debye length. 

For linear labels larger than the Debye length, McCormick et al. introduced a relation for the 
hydrodynamic drag coefficient, with which can be determined from the size of the polyelec- 
trolyte and label monomers, bpE and and the corresponding Kuhn lengths, &k,PE and &k,L> which 
describe the stiffness of the chains: 



"L = — r^^- (13) 
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Figure 4: A schematic representation of the "blob" picture used to derive Equation 13. 

The derivation of Equation 13 assumes that the polyelectrolyte and the label can be represented 
by a series of hydrodynamically equivalent entities, called "blobs" as shown in Figure Figure 4. 
The number and the size of these blobs depend on the bond length and flexibility of the chains, 
resulting in the presented relation for (Zl- 

The total effective friction of the polyelectrolyte-label compound with the surrounding solvent 
is linear in the total number of hydrodynamically equivalent monomers given by 



This is true for long polyelectrolytes in the free-draining regime, where the size of the compound 
is larger than the Debye length Xo, since the hydrodynamic interactions between the individual 
monomers are screened on this length scale, as shown in a prior study on free-solution elec- 
trophoresis. ^^'^^ 

Thus, the hydrodynamic drag (Xl can be directly calculated from the persistence lengths of the 
polyelectrolyte and of the label using Equation 13. Here, /p,PE and /p,L are calculated from the bond 
correlation function: 



= A^PE + ocNl. 



(14) 




(15) 
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where bi is the /-th bond vector and b is the average bond length. The angular brackets (...) denote 
an ensemble average.^ 

Under the chosen conditions, the persistence length of the polyelectrolyte is found to be 

/p,PE = 5.1±0.3, 

and the label's one 

/p,L=1.9±0.1. 

The difference between these two values is due to the electrostatic repulsion between the monomers 
of the polyelectrolyte. In our model, all monomers have the same size, so that Equation 13 is 
reduced to 

aL= (0.37 ±0.03) L. (16) 

'p,PE 

The comparison between the measured drag coefficient and the theoretical prediction in Fig- 
ure Figure 5 shows an agreement for the respective regimes of validity. The agreement between 
theory and experiments has been shown before. io,i8,39 p^j. ^^bels with a hydrodynamic size smaller 
than the Debye length, i.e., Rh < Ad, Equation 11 gives the correct prediction for the drag coeffi- 
cient ttL. Longer labels, however, can no longer be seen as a single polymer coil with a hydrody- 
namic size Rh, but instead the blob picture described by Equation 13 has to be used. This prediction 
is only valid when the hydrodynamic size becomes larger than the Debye length. For i?h ~ the 
expected cross-over between these regimes is observed. 

It remains to be emphasised that, by determining Ol from the measurements of the persis- 
tence lengths and the hydrodynamic radius, there is no free fitting parameter and the quantitative 
agreement in Figure Figure 5 is noteworthy. 

^For a discussion about different ways to determine the persistence length in computer simulations please refer 

to. 38 
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Figure 5: The hydrodynamic drag coefficient (Zl for the linear label is compared to the theoretical 
predictions of Equations 1 1 and 13. Both predictions are valid in different regimes with the division 
being ~ Ad indicated by the vertical dashed line. Note that neither theory has a free fitting 
parameter used to achieve the quantitative agreement with the simulation data. 
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Increasing the hydrodynamic drag coefficient 

In Figure Figure 5, it is shown that the total drag coefficient Ol for linear labels can be increased 
by using longer labels, and that beyond the Debye length the increase is linear with the length L of 
the label. Unfortunately, the experimental requirement of strict monodispersity of the label limits 
the size of linear polymeric labels that can be synthesized and prepared. In this section, it will be 
shown how one can influence the total drag coefficient also by modifying the relative stiffness of 
the polyelectrolyte and drag-tag chains. 

Increasing the label stiffness 

Equation 13 shows the dependence of ocl on the persistence lengths of polyelectrolyte and label. 
Therefore, Ol can be increased by either increasing the persistence of the label or decreasing the 
persistence of the polyelectrolyte. Both ways will be investigated in this subsection. 
First, an additional harmonic bond angle potential. 



is added to the interaction between the label monomers, where is the angle between two consec- 
utive bonds. Here, ^ba = 30 and 0o = are chosen. 

The bond angle potential increases the hydrodynamic radius of the 30 monomer label from 
i?h,L = 3.00 ±0.05 to 



and thus puts the label size into the regime where the blob picture is valid. The increased stiffness 
doubles the persistence length of the label from /p^L = 1.9 ± 0.1 to 




(17) 



i?h,L = 5.25 ±0.05, 



/p,L = 4.0±0.1, 
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which yields an increased drag coefficient according to Equation 13 of: 



aL = (0.79 ±0.04) 23.7. 
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Figure 6: The hydrodynamic drag coefficient of a stiff linear label is higher than that of a fully 
flexible label with the same length. The slowdown of the stiff label is correctly predicted by 
Equation 13 (solid line). 



Figure Figure 6 compares the theoretical predicted slowdown of the stiffer 30 monomeric la- 
bel to the measured mobilities. As before, an excellent agreement to the theory is found for the 
investigated label lengths. Please note that there is no fitting parameter here. 
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Reducing the polyelectrolyte stiffness 

In practice, it is not easy to increase the stiffness of the drag-tag in order to to increase Ol since 
this implies changing the chemistry of the label. However, the value of ttL is a relative measure 
of the stiffness of the two components of the hybrid polyelectrolyte-label molecule (see Eq. 13). 
Therefore, an increase of ttL can also be achieved if we reduce the stiffness of the polyelectrolyte. 
The persistence length of a polyelectrolyte can be reduced very effectively by increasing the ionic 
strength of the buffer: this increases the screening of the electrostatic repulsive interactions that 
stiffen the backbone of the charged polymer. In this subsection, we investigate the effect of chang- 
ing the concentration of salt from cs = 16mM to cs = IM, which reduces the Debye length from 
Ad ~ 4.2 to Ad ~ 0.65. 

The increased electrostatic screening reduces the extension of the polyelectrolyte chain and 
reduces the contribution of electrostatics to the persistence length, which is determined to be 

/p,PE = 3.8±0.2, 

whereas the label persistence is unaffected. Thus, one predicts a hydrodynamic drag coefficient of 

(XL = (0.50 ± 0.04) L, 

using Eq. 13, instead of 

aL = (0.37 ±0.03) L 

for the previous salt concentration (see Eq. 16). 

The change in electrophoretic mobilities for labels of lengths L = 20 and L = 40 can be seen 
in Figure Figure 7. Please note that the free-draining mobiUty jUq also changes due to the fact that 
the additional salt also increases the screening of the polyelectrolyte charge, thus reducing the net 
force from the external field: 

jUo = 0.1 13 ±0.002. 
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Figure 7: In the presence of 1 Mol additional salt, the persistence length of poly electrolyte is 
reduced, changing the relative hydrodynamic drag of the label. The observed mobility of the 
polyelectrolyte molecules for two linear labels of length 20 and 40 are compared to the prediction 
using Equation 13. 
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The reduction in the absolute mobility jUq together with the increase of the diffusion coefficient 
due to the more compact conformation of the molecules at higher salt concentrations negatively 
affect the size-selectivity as can be seen when evaluating the resolution factor R as defined by 
McCormick et al. If we keep the poly electrolyte length A^, the label length L, the electric field 
strength E and the elution distance constant then one obtains for a given Ol a resolution factor R 
that only depends on the diffusion coefficient Dq and the free-draining mobility jUq: 

R ~ VDo/jUq. 

The way the resolution factor is defined a higher value indicates a lower size-selectivity and thus 
the size-selectivity is decreased if the relative drag of the label is increased by increasing the salt 
concentration. Consequently, an increased hydrodynamic drag coefficient is less effective when 
achieved by adding additional salt. 

Branched drag tags 

In this section, we will investigate the use of branched labels as a possible way to synthesize more 
efficient ELFSE drag-tags. First, the results obtained in a recent experimental study by Haynes et 
al. are briefly reviewed. The study compared a linear polypeptide drag-tag with 30 repeat units 
to two branched drag-tags, each with 5 side-chains spaced evenly along a 30 unit-long backbone. 
The two different branched labels had 4 and 8 monomer long side-chains. The drag coefficients 
Ol were obtained by measuring the mobility of two different DNA fragments of 20 and 30 bases in 
length. It was found that the value of ocl increases roughly linearly with the total molecular weight 
of the branched label. 

This astonishing observation was theoretically analysed by Nedelcu et al. It was shown that 
the drag coefficient is directly related to the hydrodynamic radius (as one would expect from the 
blob picture), and that the linear dependence on molecular weight is only approximately true in the 
limit of short side chains. 
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As a matter of fact, the drag provided by a linear label is always higher than that provided by 
a branched label of the same molecular weight. The reason for this is that, with a fixed length 
backbone, a branched polymer is essentially a compact star polymer with a smaller hydrodynamic 
size than the linear equivalent. Indeed, as the number of arms increases, the branched polymer 
becomes even more compact and less favorable for ELFSE. 

Based on the observations, the following optimal design using branched polymeric labels for 
ELFSE was proposed: I) side chains with length comparable to the distance between branching 
points, or II) two long branches located near the ends of the molecule's backbone. 

Here, the focus will be on investigating the effect of the length of the side chains for a polymeric 
drag-tag with a fixed backbone length. Similar to the structure of the label used by Haynes et al. , 
the label has a backbone of L = 30 monomers to which 5 side chains are attached evenly spaced 
along the backbone. The side-chain length is varied from 2 to 8 monomers, so that the total number 
of monomers in the label ranges from 40 to 70. The drag coefficient of the labels is determined by 
measuring the electrophoretic mobility of polyelectrolyte chains from N = 20toN = 60. 

Figure Figure 8 shows the simulation results for a 30 monomeric label without side chains and 
with the tetra and octamer side chains. To analyze the hydrodynamic drag of the branched labels 
in detail, (Zl is determined according to Equation 6. The obtained values are compared to the 
corresponding value of a purely linear drag-tag with the same number of monomers. 

Figure Figure 9 confirms the work by Nedelcu, showing that the label with the highest drag 
per monomer is the linear label. For the same number of monomers L, the hydrodynamic drag 
coefficient Ol of the linear label is higher than that of the branched one. But it also shows that 
the addition of side chains can be used to increase the hydrodynamic drag of the label. This is 
attributed to two effects: firstly, the hydrodynamic size of the label is increased as the side chains 
extend from the label. Of similar importance is the second effect, namely that the side chains stiffen 
the label due to steric repulsion with the backbone, increasing the overall persistence length and 
increasing the linear length of the backbone. In fact, this is the main contribution to the increase for 
the side chains of two and four monomers as we confirmed by measuring the change in persistence 
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Figure 8: The reduced mobility /i//io for poly electrolytes with an attached linear label with five 
side chains of length 4 and 8 shows a more pronounced slowdown than for the label without side 
chains. 
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Figure 9: The hydrodynamic drag coefficient (Xl of a branched polymeric label is compared to 
the previously determined drag of a linear label. L is the total number of monomers. As long as 
the hydrodynamic radius of the label is smaller than the Debye length Ad, the (Zl is given by 
Equation 11. The vertical lines indicate the number of monomers L for which R\i{L) ^ Ad obtained 
from simulations. 
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length of the backbone. 

Interestingly, the drag coefficients obtained for the labels show a scaling with the hydrodynamic 
radius Rt, as given by Equation 11. Since the polymer coil formed by the branched label is more 
compact, it is less penetrated by ions and, therefore, the prediction of Equation 1 1 remains valid 
for a higher number of monomers compared to the linear label. 

The experimentally observed linear scaling with L can be attributed to seemingly linear rela- 
tionship between and L, but, as Nedelcu et al. have shown before, this is only approximately 
true in the case of side chains smaller or equal to the spacing along the backbone. The only relevant 
quantity in all cases is the hydrodynamic radius and its contribution to the hydrodynamic drag, as 
formulated in Equation 11. 

Although linear labels remain preferable as long as the pure hydrodynamic drag coefficient 
ttL per molecular weight is concerned, branched polymers offer practical advantages because of 
the possibility of synthesizing larger and somewhat stiffer monodisperse molecules in a simple, 
stepwise way. 

Micellar drag tags 

Recently, Grosser et al. ^^'^^ proposed another promising class of drag-tags that in principle can 
provide very large hydrodynamic drag coefficients Ol. They used nonionic surfactant Triton X- 
100 micelles that attach to PNAA-taged (PNA amphiphile) DNA strands. The micelles are water- 
soluble and are created and destroyed on a timescale of milliseconds to seconds, forming a fairly 
monodisperse populations of structures with a tunable size and morphology. During the whole 
electrophoresis time, a single DNA strand attaches to a large number of different micelles. Of 
importance for the ELFSE application is the fact that this leads to an averaging effect between 
micelles of different sizes for the individual DNA strand, meaning the DNA can be though of as 
having a drag tag of fixed size (R), where (R) is the average micelle size. Only with this averaging, 
the natural polydispersity of the micelles is overcome and a measurement with a size resolution up 
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to a couple of base pairs is possible. 

As a free DNA strand quickly attaches to a new micelle, the DNA is bound to a micelle most 
of the electrophoresis time. Consequently, the transiently bound micelles provide about the same 
hydrodynamic drag as a covalently bound drag-tag of similar size would provide. The reported ttL 
values range between 33 and 58 for a single micelle, depending on the micelle type and the PNAA 
molecule used for connecting to the DNA strand. Savard et al. showed that dual-tagging of the 
DNA, i.e., attaching a PNAA molecule to both ends of the DNA strand so that two micelles are 
transiently bound can increase the hydrodynamic drag even further. 

In this study, four different micelles with radius ^ = 2 to = 5 are attached to polyelectrolyte 
chains of different lengths. Neither the attaching and detaching process nor the forming of the 
micelles themselves are modelled explicitly; only the hydrodynamic drag of a covalently bound 
spherical drag-tag is investigated. The results by Grosser and Savard show that the polyelectrolyte 
is in fact attached to micelles for most of the time and only spends a small fraction of time without 
a drag-tag. As the micelle size can be exactly chosen in simulations, the averaging procedure 
resulting from the attaching and detaching process does not need to be included. 

Figure Figure 10 shows that micellar drag-tags can be successfully used for electrophoretic 
separation of polyelectrolyte chains. The values for ttL are obtained as before and compared to 
Equation 11, which correctly predicts the observed behaviour. With the chosen micelle radius of 
R = 5, drag coefficients up to ttL = 24.3 are achieved. 

The results show that the hydrodynamic drag is directly related to the size of the micelle, as 
can be seen in Figure Figure 11. A linear increase with the radius is observed, as expected from 
Stokes theory. Again, Equation 1 1 correctly predicts the drag coefficient, clearly indicating that 
only the hydrodynamic size of the drag-tag is important, not the number of units it is made of or 
the weight associated with it. 

The authors strongly believe that the use of micellar drag-tags has great potential for the further 
advancement of end-labeled free-solution electrophoresis. Especially the tunable size makes them 
ideal candidates, as the drag coefficient can be optimised to the lengths of polyelectrolyte frag- 
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Figure 10: The electrophoretic mobility /i as a function of the polyelectrolyte length A'^ becomes 
size dependent when a micellar drag-tag is attached. The magnitude of the slowdown depends on 
the radius R of the micelle. 
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Figure 1 1 : The effective hydrodynamic drag coefficient of a micellar drag-tag is directly pro- 
portional to its radius R. Equation 11 (the solid line) gives a very good prediction of the drag 
coefficient for all tested micelles. 
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ments to be analysed. Currently, we are investigating the usability of cationic micelles as drag-tags 
carrying positive charges, which create an additional force on the polyelectrolyte-label compound, 
possibly enhancing the size separation. 

Conclusion 

In this paper, we have presented a detailed study of end-labeled free-solution electrophoresis (us- 
ing various hydrodynamic drag-tags) by coarse-grained molecular dynamics simulations. Linear, 
branched and micellar drag-tags were investigated. The simulations support the theoretical predic- 
tions and can be matched quantitatively to it. This enables the use of computer simulation as a tool 
to support the design of improved hydrodynamic drag-tags usable for electrophoretic separation of 
polyelectrolytes in free-solution. 

It was specifically shown that the drag coefficient of the label is determined by its hydrody- 
namic size and not by its weight. The hydrodynamic drag per label monomer is largest for linear 
labels, but experimental restrictions in the synthesis of such labels and the monodispersity require- 
ment limit their practical applicabiUty. 

The addition of side chains to a linear label offers the possibility to increase the hydrodynamic 
size without having to increase the linear length of the label. The synthesis process creates perfectly 
monodisperse labels. It was shown that the label efficiency is increased with the length of the side 
chains for the drag-tag sizes studied in this work. In addition to increasing the lateral size of 
the drag-tag, the side chains also increase the persistence of the backbone and thus contribute in 
two different ways to the increased hydrodynamic size. Especially the steric stabilisation of the 
Unear backbone is responsible for an initial increase of the drag-coefficient with the total number 
of monomers of the label, i.e., with the molecular weight. For longer side chains, the lateral 
contribution to the hydrodynamic radius becomes more important. 

The third class of labels investigated seems very promising for the future of ELFSE. Transiently 
bound micelles provide a significantly higher hydrodynamic drag, as they can be prepared with 
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a large hydrodynamic radius. Additionally, the time averaging by attaching to many different 
micelles over the electrophoresis time span helps to meet the monodispersity criteria. This study 
showed that the hydrodynamic drag is directly proportional to the hydrodynamic radius of the 
micelle. The efficiency of this method is, in principle, only limited by the size of labels that can be 
prepared. 

Our results demonstrate convincingly that theory and computer models can support the experi- 
mental progress towards the design of novel improved drag-tags, thereby extending the applicabil- 
ity of the ELFSE technique. The usability of charged drag-tags is currently under investigation. 
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